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Abstract: The phase of inflationary expansion in the early universe produces superheavy 
relics in a mass window between 10 12 GeV and 10 14 GeV. Decay or annihilation of these 
superheavy relics can explain the observed ultrahigh energy cosmic rays beyond the Greisen- 
Zatsepin-Kuzmin cutoff. We emphasize that the pattern of cosmic ray arrival directions seen 
by the Pierre Auger observatory will decide between the different proposals for the origin of 
ultrahigh energy cosmic rays. 

PACS Nos.: 98.70.Sa, 98.70.-f, 95.35.+d, 14.80.-j 

Resume : French version of abstract (supplied by CJP if necessary) 
[Traduit par la redaction] 

1. Introduction 

The nature and origin of the dark matter in the universe and the origin of ultrahigh energy cosmic 
rays (UHECR) are certainly two of the most interesting problems of astroparticle physics, and maybe 
even of physics, at the beginning of the 21st century. 

The three puzzles related to the origin of cosmic rays with energies beyond the Greisen-Zatsepin- 
Kuzmin cutoff, E > E GZ k — 4 x 10 19 eV, are: 

1. Scattering off cosmic microwave background photons limits the penetration depths of charged par- 
ticles at these energies to distances < 100 Mpc [1, 2, 3]; 

2. the distribution of arrival directions of UHECRs does not seem to favor any known astrophysical 
sources within the GZK cutoff length; 

3. it seems extremely difficult to devise sufficiently efficient astrophysical acceleration mechanisms 
which could accelerate particles to energies E > Eqkz, and at the same time overcome collisional 
and radiation losses. 

The so called top-down models of UHECRs combine both problems by proposing that ultrahigh 
energy cosmic rays arise in the decay [4, 5, 6, 7, 8] or collisional annihilation [9, 10, 1 1] of superheavy 
dark matter particles. We will use both the acronym SHDM and WIMPZILLA 1 for superheavy dark 
matter particles. 

SHDM decay can arise through direct decay of relic particles ("WIMPZILLA" decay) or through 
the annihilation of superheavy relic bound states ("WIMPZILLIUM" decay), but the underlying decay 



The word WIMPZILLA was coined by Kolb, Chung and Riotto in their investigations of possible origins of superheavy dark 
matter particles [12], 
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mechanism plays no role for the predicted pattern of UHECR arrival directions. This pattern differs 
strongly from the anisotropy pattern predicted by collisional WIMPZILLA annihilation. The reason 
for the difference is that the anisotropy pattern predicted by decay models is proportional to the WIM- 
PZILLA or WIMPZILLIUM density nx (r), and is dominated by the smooth background halo. Collisional 
annihilation, on the other hand, predicts that the anisotropy pattern should be proportional to n x (r), 
and is constrained by unitarity limits on annihilation cross sections. Therefore collisional annihilation 
can only work in dense cores of dark matter substructure in the galactic halo. These dense cores were 
denoted as WIMPZILLA stars [11]. As a consequence, the collisional annihilation scenario predicts a 
pointlike source distribution with increasing density towards the galactic center. 

All modern acceleration ("bottom-up") models for UHECR origin assume powerful extragalactic 
sources, and therefore predict that the pattern of observed arrival directions should not correlate with 
the galactic halo. The different patterns of UHECR arrival patterns predicted by all the different source 
models imply that a dedicated large statistics experiment like the Pierre Auger observatory can easily 
identify the correct model from its anisotropy signal. 

Sec. 2 explains the origin of a mass window for the generation of relic superheavy dark matter 
particles during inflation. Sec. 3 summarizes and updates the calculation of the ultrahigh energy cos- 
mic ray flux from collisional WIMPZILLA annihilation, and the origin of the unique pattern of arrival 
directions predicted by that model. In Sec. 4 we compare the different anisotropy signatures predicted 
by all contemporary proposals for the origin of ultrahigh energy cosmic rays, and Sec. 5 contains our 
conclusions. 



2. Superheavy dark matter from inflation 

Particle creation as a consequence of non-adiabatic expansion was discovered already in the late 
1950s and 1960s, see [13] and references there. However, the effect was found to be negligible in 
radiation or dust dominated epochs, and therefore this mechanism for non-thermal particle creation 
was rediscovered and garnered much more interest only after the necessity for inflation was realized. 
The effect is usually considered in terms of the Bogolubov transformation between in and out vacua in 
an expanding universe [12, 14, 15, 16, 17]. Particle creation during preheating after inflation can also 
arise as a consequence of a direct coupling between the inflaton and other matter fields [18, 19]. 

Here we will expand on a simple discussion given in [11] to see how particle production in an 
inflationary universe can be understood by studying the evolution equations of weakly coupled scalar 
fields in the expanding universe. The reasoning outlined here is not intended to compete in any way 
with the traditional operator methods to study particle production from spatial expansion, but it may 
provide a complementary and helpful view. 

If interactions with other matter fields can be neglected, a scalar field in a Friedmann-Robertson- 
Walker background with metric 

d" 1 * 2 , Jj,(2 , „2 • 2„q 2\ _ jj.2 , „2n\f J„2 , /.. ( X " ^0 



ds 2 = -dt 2 + a 2 (t) ^ - _ ^ 2 + r W + r 2 sin 2 $ dip 2 j = -dt 2 + a 2 (t) [dx 2 + k K - _ ^ j (1) 
satisfies 

0(x, t) + 3^|^(x, t) - -Jj^ [(5 ij - fccV) didj - Mx'dil 0(x, t) + m 2 0(x, t) = 0, (2) 
and the corresponding energy density per unit of comoving volume is 

hOO 



e(x,t) = V^?T U0 (3) 



a 



(*) 



2V1 - kr 2 



1 



2 (x,t) +mV(x,t) + (S ij - fczV) d i( p(x,t) ■ dj4>{x,t) 
a (t) 
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The violation of time translation invariance in (1) implies violation of energy conservation, of course: 

(4) 



g(x,t) - a(t)d., 

a 2 (t)a(t) 



§v _ kx l x 3 
X ' ^ ft? 



2V1 - kr 2 



A priori the sign of the expression on the right hand side is indefinite, e.g. a rapidly evolving field of 
low mass and with small spatial fluctuations loses energy during spatial expansion as long as it evolves 
rapidly enough. However, Eq. (2) tells us that spatial fluctuations (which should never have been large 
anyway for the Friedmann-Robertson-Walker ansatz to work) are soon negligible in the inflationary 
expanding universe a(t) oc exp(iJi). This leaves us with the simple equation 

4>{t) + 3H<j>(t) + m 2 (t>{t) ~ 0. (5) 

Approximately constant H during inflation then yields for the time evolution of the comoving energy 
density 

g(t) ~ l -a%t) (j> 2 (t)+m 2 cj> 2 (t)) 

~ A+ cxp(tV9H 2 -4m 2 ) + A_ cxp(-ty/9H 2 - 4m 2 ) + B. 

This implies a growing mode in the comoving energy density of weakly coupled states with m < 
1.5H ~ 10 14 GeV. What is special about the superheavy particles is that their comoving energy density 
is conserved after inflation, because the behavior of massive (m > t^ 1 ) weakly coupled states in the 
subsequent radiation and dust dominated backgrounds preserves their energy. The asymptotic solution 
for weakly coupled massive states with m > i -1 in such a background yields (with i = 3 for dust, 
£ = 4 for radiation) 

4>(t) oc t~ 3/e cos(mt + ip) , 
g(t) oc a 3 (t)t-^ e oc t°, 

and this implies in particular that the comoving density of massive particles freezes out at the end of 
inflation (t ~ 10" 36 s) if 



rn 



> r 1 ~ 10 12 GeV. 



These considerations indicate a mass window for direct gravitational production of superheavy relic 
particles during inflation 

10 12 GeV < m < 10 14 GeV. 

After inflation intrinsically unstable matter of mass M will decay on time scales r < Mp [anck /M 3 , 
where the upper bound assumes that the particles couple to their decay products at least with gravita- 
tional strength. This means that intrinsically unstable superheavy particles should not be relic particles. 
This difficulty has motivated the collisional annihilation scenario for ultrahigh energy cosmic rays from 
superheavy dark matter [9], because the unitarity limit on reaction cross sections 



M < — (6) 



47T 

implies that superheavy dark matter without direct decay channels will still be around. 
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3. Calculation of the flux 

We consider decay or annihilation of superheavy dark matter particles of mass M x > 10 12 GeV. 
The spectral fluxes at our location r© from decay or collisional annihilation of the dark matter particles 
of density nx (r) are then 

dM(E,M x ) 1 n x (r) (?) 



47r|r - r| 2 r d 



and 



tW(E,2M x ) f 3 v 2 

Ja{E) = M J d r i 67 r|r -rP ^ (r) M > 



(8) 



respectively. Here dN (E, Ei n ) is the number of particles in the energy interval [E, E + dE] emerging 
from a decay or annihilation event of initial energy Ei n . dN{E, Ei n )/dE is related to fragmentation 
functions via 



<W(E,E in ) _ sr 1 Ar« J_ VpWr7 . 



d£ ^ CTA cLE 

where a; = E/Ei n and (a;, S 2 ^) is the differential number of particles of species i generated in the 
prescribed x-range in a decay or annihilation event with s = Ef n . The factor v in Eq. (8) equals 4 if 
the X-particles are Majorana particles, and 1 otherwise. 

iV-body simulations predict that usually about / c / <~ 5 — 10% of dark matter halos should exist in 
substructure. Therefore in the decay models the ultrahigh energy cosmic ray flux will be dominated by 
the smooth background halo. Evaluating (7) e.g. for a Navarro-Frenk-White halo [20] 

nx(r) = 4:nx{r s )r 3 s /r(r + r s ) 2 

yields 

The scale radius r s for the Milky Way is not well known, but for the whole range of say 5 kpc < r s < 
50kpc comparison between (9) and the ultrahigh energy cosmic ray flux from AGASA [21] yields 
that decaying superheavy dark matter should only make a small contribution to the galactic dark matter 
density. Of course, the problem is to explain a lifetime Td > 10 17 s for unstable superheavy dark matter. 

Eq. (8) can also be evaluated analytically for a Navarro-Frenk-White halo, but the unitarity limit 
for s-wave annihilation 

(w)=l x^= S x440x 1 0- m »/sx(~) 2 xl^ ! , £<!, „0, 
M x v \ Mx J v 

implies that any ultrahigh energy cosmic ray flux from collisional annihilation in the background halo 
is negligible. Therefore collisional WIMPZILLA annihilation can only make a noticeable contribution 
to the ultrahigh energy cosmic ray flux if it originates in dense cores of dark matter subclumps of 
the galactic halo [9]. These dense cores were denoted as WIMPZILLA stars in [11], and a simplified 
estimate of the cosmic ray flux from WIMPZILLA stars yields 

. /m N cl V core dAf(E,2M x ) 2 0.1f cl M halo dAf{E,2M x ) . . 
3a{E) " V ^dP dE nxM * V MM X dE Hx{aAV) - (U) 
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Here dr 2 = (r~ 2 ) is a mean inverse distance squared for our separation from galactic WIMPZILLA 
stars. In our estimate in [11] we had used d ~ lOkpc, which seems reasonable: e.g. our inverse mean 
distance squared to visible galactic substructures, the globular clusters, corresponds to d = 7.3 kpc. 

We have parametrized the unknown annihilation cross section already in terms of the unitarity 
bound (10). Now we would like to compare e.g. with the UHE cosmic ray spectrum in Ref. [9], Fig. 
1 (note that the cross section stated there would only be required for collisional annihilation in the 
smooth background halo and does not apply to annihilation in dark matter subclump cores). For the 
comparison we parametrize the core density nx of WIMPZILLA stars in terms of the solar density 



-TrrriT-v " im x 11,17,11 



This yields a flux at 10 11 GeV of order 

E*j(E) = 5.74x 10- eV 2 m" 2 s" 1 x J± x - x ( 7 -*£°) \ VX *<12) 

v ; b=ioii Gcv 0.05 2 x 10 12 M Q \ d J lO" 3 

This complies with the spectral fit for the flux per solid angle in Ref. [9], Fig. 1, if t?£ < 1, as would 
be expected for consistency of the model: If the WIMPZILLAS are not Majorana particles (y = 1) 
the required value for the fit is ri£ ~ 9.8 x 10~ 4 , and otherwise the required value would be ri£ ~ 

2.5 x 10~ 4 . 



4. Expected anisotropy patterns beyond 20EeV 

Fits of fragmentation functions to the ultrahigh energy cosmic ray spectrum indicate that in the 
collisional annihilation scenario cosmic rays above 20 EeV should be dominated by the fragmentation 
products of WIMPZILLA annihilation [9]. Therefore the expected anisotropy patterns discussed in this 
section only apply in this energy range. 

The purpose of this section is to emphasize that the anisotropy pattern observed by the Pierre Auger 
observatory above 20 EeV will provide a crucial direct test of the different proposals for the origin of 
ultrahigh energy cosmic rays: Bottom-up acceleration of charged paricles in AGNs (see e.g. [22]) or 
gamma ray bursts [23], ultrahigh energy neutrinos travelling over cosmological distances and creating 
Z-bursts [24], WIMPZILLA or WIMPZILLIUM decay [4, 5, 6, 7, 8], or collisional annihilation in WIM- 
PZILLA stars [9, 11]. 



Origin of UHECRs 


Expected anisotropy pattern 


Collisional annihilation in WIMPZILLA stars 


About 1000 pointlike sources with increasing 
density towards the galactic center. 
No correlation with galactic SNRs. 


WIMPZILLA or WIMPZILLIUM decay 


Dominated by uniform increase towards galactic center. 


Z-bursts 


Approximately isotropic distribution with only weak 
correlation to structure within 150 Mpc. 


Bottom-up acceleration 


Correlation with local superstructure within 150 Mpc. 
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Once fully operational, the Pierre Auger observatory should see more than 300 events per year with 
energies above 40 EeV. If the collisional annihilation scenario is correct, the observatory should see a 
large number of multiplets within its angular resolution, with increasing density towards the galactic 
center. 

5. Conclusions 

Collisional annihilation of superheavy dark matter particles as a source model for ultrahigh energy 
cosmic rays has the advantage of naturally explaining the absence of a GZK cutoff in the spectrum 
without correlation to local AGNs, without the need of postulating an extremely powerful and efficient 
acceleration mechanism, and without the need to explain extremely long lifetimes of intrinsically un- 
stable particles. It has the disadvantage of having to postulate formation of a few relatively dense cores 
of dark matter subhalos or subhalo remnants. 

The two different top-down scenarios of decay or collisional annihilation imply qualitatively differ- 
ent anisotropy signals in cosmic ray arrival directions, and these also differ from the anisotropy signals 
predicted by the various bottom-up scenarios. What constitutes a statistically significant dataset to con- 
firm or reject the different anisotropy patterns depends on the different patterns, of course. However, for 
the collisional annihilation scenario the Pierre Auger observatory should see a large number of multi- 
plets within its angular resolution already after about two years of full operation, because the UHECRs 
should arise from a limited number of pointlike sources in our galactic halo. Based on current parti- 
cle physics extrapolations for chemical composition of cosmic rays to very high energies, bottom-up 
acceleration scenarios seem much more popular in the cosmic ray community. However, collisional 
WIMPZILLA annihilation is not an overly speculative theory, and based on its expected anisotropy 
pattern it will be confirmed or rejected in a very short time frame. It would therefore seem prudent to 
await what the Pierre Auger observatory will tell us about patterns of UHE cosmic ray arrival directions 
before jumping to any conclusions about the validity of different models. 
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